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A new defect-assisted mechanism of multiphonon intraband carrier relaxation in semiconductor 
quantum dots, where the carrier is found in a coherent superposition of the initial, final, and defect 
states, is proposed. It is shown that this mechanism is capable of explaining the observed trends in 
temperature dependences of the intraband relaxation rates for PbSe and CdSe colloidal nanocrystal 
quantum dots. 
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I. INTRODUCTION 

Recent experimental studies of PbSe nanocrystals 
(NCs) have revealed an unexpectedly fast intraband re- 
laxation in semiconductor quantum dots (QDs) with en- 
ergy separation between adjacent electron (or hole) levels 
exceeding by far the optical phonon energies A pro- 
nounced temperature dependence of the observed relax- 
ation rates [3 suggests that underlying mechanism should 
involve multiphonon transitions while Augcr-likc relax- 
ation processes [H are ruled out due to a nearly symmet- 
ric energy spectra in the conduction and valence bands 
of PbSe. 

The multiphonon transitions between otherwise or- 
thogonal quantized electronic states in a QD may oc- 
cur due to nonadiabaticity of the electron-phonon sys- 
tem I3I1HQ- A model study performed in conjunction 
with the experiments on PbSe NCs revealed that, in order 
for this mechanism to be responsible for the experimen- 
tally observed temperature dependence of the relaxation 
rate, the electron-phonon coupling should be stronger 
than one would expect 0- However, the study of detailed 
models of electron-phonon interaction in quantum dots 
(upon which the judgement about the coupling strength 
is based) has always been a weakpoint in approaches to 
energy relaxation |8| and Raman |9| processes. In case of 
PbSe NCs the situation is further obscured by the fact 
that the widely accepted model of their electronic struc- 
ture [13 which would be natural to pick for an estimate 
of electron-phonon coupling strength fails to describe ab- 
sorption spectra of PbSe NCs[l| and is therefore not quite 
reliable. Thus, instead of seeking to either confirm or 
rule out this mechanism we will concentrate on examin- 
ing other possibilities for the fast energy relaxation. 

As it has been proposed for epitaxially grown QDs [ill, 
such possibilities can be provided by a localized state of 
an impurity or defect close to the QD surface. In case of 
the experiment of Ref. [l] such states can correspond to 
deep impurities in the silica host surrounding PbSe NCs 
occurring in an appreciable proximity to the NC surface. 
Otherwise they can be due to surface states. An electron 
trapped to the impurity is strongly localized which as- 
sures its strong coupling to local lattice vibrations. One 



can then consider a relaxation process involving the inte- 
rior QD levels as the initial (2) and final (0) states, and 
the state at the impurity as an intermediate (1) state as 
shown in Fig. 1. Although processes of this type have 
been theoretically investigated I MTllfl2l . all these studies 
considered the relaxation as a two step process, so the 
relaxation time r = T2^o = ^2^1 + ^l^o- But if ti^o 
is shorter than the phonon lifetime in state 1 then the 
transition 2 — > must be treated as a coherent quan- 
tum mechanical process. Such process represents a new 
mechanism of phonon-assisted carrier relaxation in QDs 
and its studying is the scope of the present work. Similar 
processes have been considered for capture or generation 
of electron-hole pairs by impurities in bulk semiconduc- 
tors GJGJ 

Another system where the proposed mechanism can 
play the key role is provided by CdSe colloidal nanocrys- 
tals. In the recent experiments on such quantum dots 
the fast intraband relaxation of electrons has been ob- 
served in the absence of the holes The electrons were 
injected into the quantum dots instead of being optically 
generated, and the p ump -probe measurements were sub- 
sequently carried out I la. These experiments showed that 
intraband relaxation times vary significantly with surface 
modification suggesting involvement of a surface state 
such as the empty s - state of the surface Cd 2+ atoms [H. 
A faster relaxation rate was observed for lower temper- 
atures [TH though more thorough studies of the tempera- 
ture dependences of the relaxation rate in CdSe colloidal 
nanocrystals are desirable. 

It has been noticed that multiphonon relaxation pro- 
cesses are analogous to tunneling problems in electron 
transport utilizing the tunneling hamiltonian. It turns 
out that the relaxation mechanism proposed in this paper 
is analogous to the problem of electron tunneling through 
a potential barrier placed between two leads and contain- 
ing localized electronic states coupled to phonons and 
resonant with the energy of the tunneling electron. The 
latter problem was studied by Glazman and Shekhterflfjl. 
There is, however, a substantial difference between the 
two problems due to the descrete energy spectrum for 
carriers confined in a QD. Indeed, in the problem of elec- 
tron tunneling there is a continuum of both initial and 
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final states in the leads. Thus, this problem has a char- 
acteristic time related to the elastic tunneling of the elec- 
tron from the localized state to either of the leads and 
determining the width of the resonance. The rate of the 
phonon-assisted tunneling can be then calculated using 
some generalization of the Fermi golden rule with the 
(5-function accounting for the energy conservation being 
integrated out while summing over the continuum of fi- 
nal electronic states. In our problem the initial and final 
electronic states are descrete. However, the spectrum 
of the entire system "electron plus phonons" is continu- 
ous and the scheme of calculation of the transition rate 
proposed in Ref. Qjl can still be applied, although some 
changes should be introduced. First, the width of the 
resonance will be determined by the lifetime of phonons 
in the intermediate state. Second, the phonon density of 
states will enter the expression for the transition rate. It 
turns out that these changes can be more naturally in- 
corporated using a more traditional approach based on 
calculation of overlap integrals of eigenfunctions belong- 
ing to shifted parabolic potentials (see e.g. [H). Another 
advantage of this method is that it allows one to make 
some qualitative conclusions without having done actual 
calculations. 

In this paper the proposed mechanism of carrier re- 
laxation will be studied in detail for a model QD system 
while discussions of its relevance to particular realizations 
of QDs are left for future studies. The rest of the paper 
is organized as follows. In Sec. [II] we formulate the prob- 
lem in a form most convenient for comparison with the 
work by Glazman and Shekhter [l6| and with the case of 
nonadiabaticity-induced phonon-assisted transitions be- 
tween intrinsic QD states @- The relaxation rate is ex- 
pressed through a correlation function which should be 
calculated taking into account phonon decay. In Sec. IIIII 
we devise an approximation allowing one to avoid eval- 
uation of this correlation function. The lifetime of the 
intermediate state in the relaxation process is derived in 
a consistent way. In Sec. [TV] results of Sec. OH] are imple- 
mented into calculation of the relaxation rate using di- 
rect evaluation of the overlap integrals. Dependences of 
the relaxation rate on temperature and electron-phonon 
coupling strength arc discussed in Sec. Ivl 



and 



V = £00*0 + [ei + A(fe f + b)]c ] c + e 2 Sd+ 
(Soio+c + H.c.) + (A 12 cU + H.c.) . 



(1) 



Here a, c, d are the electronic annihilation operators 
of the electronic states with the energies eq, £i, £2, re- 
spectively, and b are the phonon annihilation operators. 
We will first neglect the phonon decay and introduce it 
later on. Parameter A determines the shift between the 
minima of the vibrational potentials (corresponding to 
local vibrations and associated with different electronic 
states) whose dependence on the configuration coordi- 
nate is sketched in Fig. 1. A 12 and B i describe electronic 
tunneling between the QD states and impurity level (see 
Ref. [HI for details). 

We write the electronic wave function as 



*(t) = [u(i)c f + v{t)d) + w(t)a f ] 



(2) 



where |0) is the electronic vacuum and u(t), v(t), w(t) are 
the operators acting on the phonon subsystem. The wave 
function j2]) satisfies the Schrodinger equation in the in- 
teraction representation with respect to the Hamiltonian 
Ho 



(3) 



with initial conditions 

u{t = 0) = w(t = 0) = 0, v(t = 0) = 1. (4) 

Substituting Eqs. ©-© into Eq. (J3J), wc obtain the fol- 
lowing system of equations 

%— = [ei + V(t)]u{t) + B w w(t) + A 12 v(t) , (5) 
ot 



dw „ . . „ „ . . 
1— = £ w{t) + B 01 u{t) 
ot 



dv . . 

i-r^ = S2v{t) + A 21 u{t) , 



(6) 
(7) 



where 



II. MODEL 



The Hamiltonian describing multiphonon nonradiativc 
transitions between three electronic levels with unper- 
turbed energies eo, £1, £2 can be represented as H = 
Ho + V, where Ho is the phonon Hamiltonian describing 
local vibrations in the electronic state with the energy 
£2, 



H = aj b 1 b, 



V(t) = A(&t e luJot +be~ luJot ) . 

We will follow Rcfs. 0IH and introduce the auxiliary op- 
erator uo (t) as the solution of the following Cauchy prob- 
lem: 



fio(t = 0) = l. 



In an explicit form it reads 
A 

.LU 

A 

w 



60 = exp( — \b(e~ luJot - 1) - H.c: 

< LUo 



A Y . . o 

— (1 uot — i sin Loot) > . 

UJo 1 > 



With the help of this operator Eq. ([5]) can be written as 
d 

i -q («o u) = e\ Uq 1 u + Uq 1 [B 10 w(t) + A 12 v(t)] . (8) 

Using Eqs. the operator w(t) can be calculated 

in the lowest order in the tunneling matrix elements Bqi, 
A 12 as 
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w(t) = -Box A 12 J dt J dhe^-^e^-^x 



The probability of the electron transition is given by 



(9) 



W 2 _ = 



j_ = lim mw® 

T 2 ^0 '-* 00 t 



The average here is taken over the initial states of the 
phonon subsystem at t = 0. It is assumed that these 
states are equilibrium states. We also have taken into 
account that, at t = 0, there is no interaction-induced 
re-normalization of the phonon states since the impu- 
rity electronic level is empty. Substituting Eq. (J5]) into 
Eq. (JTUJ) we obtain 



-i{t' Q — *o)eo • 



W 2 ^ Q = lim Z/t J dt J dti J dt' J dt[ 



e <(ti-t„-ti+*o)=i e M*i-*i)0'(to,*i,«{,,«i) , (11) 

where we have introduced the correlation function 
U(t ,h,t' ,t[) = (Mt'i^it'Mto^ih)) , (12) 
and 

Z = \B 01 \ 2 \A 12 \ 2 . 

The correlation function (fT2"|) is analogous to the one 
which appeared in the problem studied by Glazman and 
Shckhter [lj|. It was evaluated in Ref. neglecting 
phonon decay. In our case phonon decay should be taken 
into account as it determines the width of the resonance 
corresponding to the intermediate state in the relaxation 
problem under study. However, the problem of evaluat- 
ing the correlation function (fT2)) taking the phonon de- 
cay into account is very difficult. Therefore, in the next 
section we devise an approximation which allows one to 
account for phonon decay without evaluating correlation 
function JT! 



III. THE WIDTH OF THE RESONANCE 

In order to understand the origin of the width of the 
resonance corresponding to the intermediate state in our 




FIG. 1: Energy diagram representing adiabatic potentials for 
the quantum dot states (0 and 2) and the impurity level (1). 

relaxation problem we will consider an auxiliary prob- 
lem which will help us to apply the concept of quasi- 
stationary states to the states of the system "electron 
localized at the impurity plus phonons" . More precisely, 
let us set Bio = A\ 2 = and consider the electron local- 
ized at the impurity and coupled to phonons. One can 
introduce the retarded Green's function of the electron 
as 

G R (t) = -i «0|c(t)ct(0) + ct(0)c(t)|0» 9(t) , (13) 

where c(t) is the electron annihilation operator in the 
Heisenberg representation and ((0|...|0)) denotes a quan- 
tum mechanical average over the electronic vacuum and 
a statistical average over the phonon degrees of freedom 
described by the Hamiltonian H$. As the electron cre- 
ation operator c^(0) enters the Green's function at t = 0, 
the latter reflects what happens to the electron when it 
is suddenly brought to the impurity site from a remote 
location (where it does not experience any coupling or en- 
tangling to phonons) at t = 0. Therefore, this auxiliary 
problem is essentially equivalent to the initial relaxation 
problem. The only difference is that, as Ai 2 = Bqi = 0, 
the auxiliary problem does not describe dynamics of the 
electron tunneling to the impurity site. The Green's func- 
tion (|13p can be written in terms of the operator mq E3 



G R (t) 



-iei t 



The spectral function denned as 1 171 

CO 

A{u>) =-23 / dtG R {t)e iut 



(14) 



(15) 



consists of a series of (5-functional peaks separated by the 
phonon energy loq. Now, if there are processes leading 
to the phonon decay then the expression for (iio(t)) be- 
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comes [18| 

(uo(*)) =cxp|-A 2 
h 



n + 1 



,-(7(wo)+«^o)< 



(7( w o) - «W ) 



(7(^0) + i^o) 2 
^ ^g-(7(^o)-«wo)< _ _|_ 

7(0) (2n + 1) - i wo 



[ 7 (0)]5 



(16) 



Here no = n(wo) stays for the Planckian factor, 7 (w) 
is the frequency-dependent phonon decay rate defined 
through the imaginary part of the phonon polarization 
operator [l8|. At zero frequency 7(0) = (which reflects 
the fact that the zero phonon absorption linewidth is not 
affected by the phonon decay) [l8lfl9l . Thus, when (p~6|) is 
substituted into Eqs. (fl4|) . (fT5|) . each peak of the spectral 
function, except for the one corresponding to no phonons, 
becomes broadened. It is natural to expect that it is this 
kind of broadening which causes the finite width of the 
resonance associated with the impurity in our relaxation 
problem. 

Assuming y(u>o) = 7 « uo, the shape of each peak 
of the spectral function can be described analytically. 
Eq. mi) yields 

(uo{t)) = exp {5 [e _7 *[(2n + 1) cosui n t — i sinwoi)] + 

iu t - (2n + 1)]} , (17) 

where the Huang-Rhys factor S is given by 
S ( w o) 

Substituting Eqs. ([T7]). (fH]l into Eq. (fTS]) we obtain 

OC 

A(uj) = ^ Ai(uj - wi), uji = 61- Sujq + luo , 

i— — 00 

Ai(u) = 2e lfiuJ0,2 ' {2Tl0+1)s x 

OO 

SR / dtIi[2S^/fi Q {n + l)cxp(- 7 t)]e Mt . (18) 



Eq. (|18p contains the averaged number of phonons, no. It 
is convenient to write out explicitly contributions of the 
initial states with given numbers of phonons. To that 
end we will use 



2/.- 



1 



fc=0 



k\{k + l)\ 



and introduce z = exp (— (3ujq) so that no = z/(l — z). 
Substituting this into Eq. (fT8"|) we obtain 



A l {io) = 2e- s Y J 



7(2fc + 



fc=0 



uj 2 + (2fc + 2 7 2 ' 



(1-z) 



2/,'+/ 



exp 



2S'z 



?2k+l 



z) k\(k + l)\ 



(19) 



This expression can be recast using the gene rating func- 
tion for associated Laguerre polynomials 120 



2S2 



6XP (-f=7) =(l-^ fc E^ L n fe+i ( 2 ^ 



(1 - 2 ) 2fc +* 



As a result, Eq. (fT9)) can be recast in the form 

00 

^(w) = (1 - e-*"°) £ A im ( W )e- m ^° , (20) 

m=0 

where 



A ; ,„(cj) = 2e s x 

^ (2A + I)7 



2k+l 



L^ l k {2S). (21) 



^ w 2 + (2/j + 2 7 2 fc!(fc + ! 



Each term in the r.h.s. of Eq. (f2"D"]) corresponds to a 
transition from the initial state with m phonons into the 
impurity state with m + I phonons and describes this 
transition's contribution into the Z-th peak of the spectral 
function Ai (oj) . The latter is contributed by all processes 
where exactly I phonons are emitted. Each of the con- 
tributions Aim(uj) consists of (m + 1) Lorcntzians. The 
Lorentzians with small k in the r.h.s. of Eq. (j2"Tj) can 
have negative weights while the weight of the Lorentzian 
with k = m is always positive. The latter Lorentzian 
determins spectral wings of Ai m (uj) and it is natural to 
associate its width, 7 (2m + I), with the linewidth of the 
transition between the states with m and m + l phonons. 
From the theory of spectral linewidths it is well known l2ll 
that the spectral width of a line corresponding to a tran- 
sition between two quasistationary states is given by the 
sum of their decay rates. As 7 is the phonon decay rate 
and we have m phonons in the initial state and (m + I) 
phonons in the impurity state, we finally arrive at the as- 
signment of the width 7 (m+Z) to the state of the electron 
at the impurity entangled to (m + I) phonons. 



IV. THE TRANSITION RATE 

We have seen that the lifetime of the intermediate state 
in the relaxation problem under study depends on the 
number of phonons associated with it. In this case it is 
more convenient to consider the relaxation process in the 
framework of an approach based on calculation of over- 
lap integrals between eigenfunctions belonging to shifted 
parabolic potentials. These potentials for our system are 
schematically depicted in Fig. [T] For the quantum dot 
states the potentials are given by@ 



Vb(<?) 



V 2 (q) 
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Mujlq 2 



£2 



(22) 
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FIG. 2: (a) The dependence of the matrix elements \M m \ on the Huang- Rhys factor, S, calculated for different values of 
m = 0,1,2,3 with I = 1, puio = 200 meV, p = 13, 7 = 0.2ps _1 . The vertical dashed line indicates the maximum of the 
dependence |Mo| 2 (5) at S = 5.7. (b) Adiabatic curves, Eqs. (|22ll - (|24fl . calculated for S = 5.7. The horizontal lines indicate 
positions of the energy levels in the potential V2 with m = and m = 1. 



Here cj is the frequency of the local vibrational mode, q 
is the corresponding configurational coordinate, and M 
is the mass of the impurity's ion. The adiabatic potential 
for the intermediate state is given by@ 



E 



(m - n)oj + e 2 - £1 + Suj + ijn 



(28) 



Vl(q) 



where 



'in 



= y/2S/{Mu) , 



■ei, 



(24) 



(25) 



£1 = £1 - S w 



The energy (5-function in Eq. (|26|) has been integrated out 
yielding the phonon density of states o;^ 1 while the value 
of uj Q is adjusted to make the number p = (e 2 — £o)/ w o 
integer. We also took into account that, according to 
Sec. IIII1 the width of the intermediate state is given by 
the number of phonons in that state times the phonon 
decay rate. 

The overlap integrals entering Eq. ([28]) can be ex- 
pressed in terms of the generalized Laguerre polynomi- 
als [111 



Let us introduce the eigenfunctions ipn\q), ipn'(q), 
and ip. 



(0) 



(q) in the adiabatic potentials V%(q), Vi(q), 



and Vo(q), respectively. These are eigenfunctions of 
the harmonic oscillator satisfying the following relation 



-S/2, 



min(m,n)! /gl m -™l j\m- 



ipn\q) = 4>n'{q — qo) = ipn'(q). Then the transition 
rate is given by the Fermi golden rule 



',(1) 



(o), 



max(m, n)\ 



i(m,n) 



(S). 



w 2 ^ a = 



2ttZ 
u> 



2—0 



where 



lM a 



(1 



-f3uj 



) J2 e- 0mMo \Mr, 



(26) 



(27) 



The relaxation is most efficient when the energy of the 
intermediate state is in resonance with that of the ini- 
tial and final states. We will assume that the resonance 
condition is always fulfilled and, therefore, 

£2 — £1 + Sujq = Iujq, leN. 

In this case the term with n — I + m becomes dominant 
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in Eq. (|28|) , and we obtain 

Mm « -tt^-t^I^SI)^!^) ■ (29) 



V. RESULTS AND DISCUSSION 

At zero temperature, only the contribution with m = 
survives in Eq. (|12T|) , and the dependence of the relaxation 
rate on the strength of the electron-phonon coupling is 
governed by the dependence of |Mo| 2 on the Huang-Rhys 
factor S. This dependence is shown in Fig. [21 (a) (solid 
curve). Both S and e\ were changed in such a way that 
the energy difference Ei — e\ + Su>q = Iloq remained fixed. 
This is equivalent to a horizontal shift of the adiabatic 
curve 1 in Fig. [21 (b). As one can see from Fig. [21 (a), the 
dependence |Mo| 2 (S') has a maximum yielding the opti- 
mal value of the Huang-Rhys factor at which the relax- 
ation is the most efficient. This can be understood using 
quasiclassical language. According to the quasiclassical 
theory, the overlap integral of cigenfunctions belonging 
to different adiabatic potentials is maximum at those en- 
ergy values where the adiabatic curves intersect @. The 
optimal value of the Huang-Rhys factor corresponds to 
the situation where both intersection points are close to 
the energy level with m = 0, as shown in Fig. [21 (b). This 
figure also explains why maximum of the curve |Mi| (5), 
also shown in Fig. [21 (a), is close to the maximum of 
|Mq| 2 (iS). As temperature is increased, thermal activa- 
tion leads to increasing importance of the terms in the 
r.h.s. of Eq. (|2"T|) with m > 0. In the meantime, relative 
contribution of the term with in = to the transition rate 
decreases due to the factor of 1 — exp(— f3uj ) in Eq. (|27|) . 
Therefore, the character of the temperature dependence 
of the relaxation rate is determined by the ratio of \Mq\ 
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FIG. 3: Temperature dependences of the relaxation time 
calculated for different values of the impurity level position 
{I = 1,2,3,5,7) with Z = (2.5) 4 meV 4 , S = 5.7. Other pa- 
rameters are the same as in Fig. [2] 



and \Mi\. One can see from Fig. [21 (a) that, for / = 1 
and Huang- Rhys factor close to its optimal value, \M\\ 
exceeds |Mo|, and relaxation time decreases with tem- 
perature as shown in Fig. 3 (solid curve). For I > 1 in- 
tersection points would be closer to the level with m = 
than to the level with m = 1 resulting in \Mq\ > \M\\. 
This will lead to increasing temperature dependence of 
the relaxation time, as shown in Fig. 3. 

Therefore, the present mechanism can reproduce both 
increasing and decreasing temperature dependences of 
the relaxation time. Those are two types of temperature 
dependences observed experimentally Note that 

when temperature sweeps from zero to 300 K, relaxation 
time in Fig. 3 changes within one order of magnitude 
which is consistent with experimental observations I2IIT5I. 



VI. CONCLUSIONS 

We have proposed a new defect-assisted mechanism of 
multiphonon intraband carrier relaxation in semiconduc- 
tor QDs, where the carrier is found in a coherent su- 
perposition of the initial, final, and defect states. Com- 
bination of different approaches has helped us to de- 
vise a controlable approximation for the relaxation rate. 
We have shown that the proposed mechanism is capable 
of reproducing both decreasing, as it was observed for 
PbSe NCs0, and increasing, as it was reported for CdSe 
QDsll5l. temperature dependences of the relaxation time. 
For reasonable values of parameters, the change in the 
relaxation time with temperature rise from cryogenic to 
room temperatures was found within one order of magni- 
tude, in agreement with experimental observations I211 1 a. 
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